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INTRODUCTION 


In  recent  years  the  Army  has  devoted  considerable  attention  to  the  use  of 
pneumatic  devices  (airbags)  in  reducing  the  landing  shock  of  various 
air-delivery  systems.  A  number  of  different  designs  have  been  investigated 
experimentally  by  Nykvist1,  Patterson2  and  others.  The  fundamental  theory  was 
presented  first  by  Browning3,  then  extended  by  Esgar  and  Morgan 4 . 

The  present  report  describes  a  preliminary  study  of  possible  systems  for 
automatic  control  in  the  course  of  a  platform  landing  cushioned  by  a  simple, 
cylindrical  airbag.  The  original  impetus  for  the  study  was  a  discussion 
between  Dr.  C.K.  Lee  and  the  author,  in  which  Dr.  Lee  mentioned  that  small 
motors  are  now  available  that  have  response  times  in  the  millisecond  range. 

The  possibility  of  using  such  a  motor  to  drive  a  device  for  opening  and 
closing  the  vent  of  an  airbag  led  to  the  investigation  reported  here. 

The  basic  theory  of  airbags  is  reviewed  in  the  next  section  and  put  into  a 
dimensionless  form  slightly  different  from  (but  wholly  equivalent  to)  that  of 
Browning3.  The  resulting  nonlinear  system  of  three  ordinary  differential 
equations  cannot  in  general  be  integrated  exactly.  However,  Section  3 
describes  two  solutions  in  closed  form  that  can  be  obtained  under  certain 
assumptions  about  the  vent  area.  Several  other  plausible  control  systems  are 
discussed  in  Section  4,  and  Section  5  presents  the  results  of  numerical 
solutions,  comparing  the  performances  of  these  control  systems.  The  results 
are  discussed  and  conclusions  presented  in  Sections  6  and  7. 


2.  BASIC  AIRBAG  EQUATIONS 


This  section  is  closely  related  to  the  analyses  given  by  Browning3  and 
Esgar  and  Morgan1*.  The  context  and  basic  assumptions  are  described  below. 

a.  The  load,  W,  is  descending  steadily  beneath  a  parachute  and  is 
attached  to  the  top  of  a  plane  platform,  dropping  vertically  and  oriented 
horizontally,  during  the  entire  impact  process.  The  initial  velocity  of 
descent  is  Vo. 

b.  An  airbag  in  the  form  of  a  cylinder  with  a  vertical  axis  and 
horizontal  end  surfaces  is  attached  to  the  underside  of  the  platform.  The 
airbag  has  height  H  and  cross-sectional  area  AB  and  initially  contains  air  at 
atmospheric  pressure,  pa.  There  is  a  vent  to  the  atmosphere  with  area  Ay. 

c.  At  time  t  =  0  the  lower  face  of  the  bag  makes  contact  with  the  ground. 
Thereafter,  the  bag  preserves  always  the  same  cross-section  area,  Afi,  but  the 
volume  decreases  as  the  platform  descends. 


d.  The  pressure,  p,  and  mass-density  p,  of  the  air  in  the  bag  are  related 
by  the  ideal  gas  law, 

p/pa  =  (e/pa)Y  (1) 

where  pa  is  the  density  of  air  at  atmospheric  pressure  and  y  is  the  ratio  of 
the  specific  heats,  specifically  y  =  1.4  for  air. 

With  these  assumptions  the  equations  of  motion  can  be  written  as  a  system 
of  three  nonlinear  ordinary  differential  equations  in 
y  =  height  of  platform  above  ground 
V  =  velocity  (positive  upward) 
p  =  mass  density  of  air  in  bag 

namely 

dy/dt  =  V  (2) 

dV/dt  =  g  [-1  +  <D+Rb)/WJ  (3) 

Abd(py)/dt  =  -CvAvpQq  (4) 


-2- 


where 


g 

acceleration  of  gravity 

D 

=  canopy  drag 

rb 

=  force  transmitted  to  platform  from  airbag 

q 

-  speed  of  air  flow  from  vent 

Cv 

=  vent  flow  coefficient 

further 

D 

=  iPaAcv2cD 

Ac 

=  drag  area  of  canopy 

CD 

=  drag  coefficient  of  canopy 

rb 

=  (P-Pa)AB 

q2 

=  JoV2  +  S 

(5) 

(6) 

(7) 

(7) 

(8) 


Equation  (8)  is  Bernoulli's  law  for  the  flow  of  air  out  the  vent  or 
orifice.  Jo  is  a  constant  that  affects  the  definition  of  the  flow  upstream  of 
the  vent  in  Bernoulli's  'aw,  and  S  has  a  form  which  depends  on  the  pressure  in 
the  bag.  Let 

Pc  =  critical  pressure  =  pg[ 1  +  ( y 1 ) /2  1 ^  =  1.893  pg  (9) 

Then 

S  =  [2 y/(Y-l))  [(p/p)  -  (P a/Pa>]  if  P  <  Pc  (10) 

=  fiPa/pa)(P/'Pc)Y/,(Y'1)  if  P  >  PC  (H) 

Equation  (2)  is  merely  the  definition  of  velocity,  conservation  of 
platform  momentum  is  embodied  in  Equation  (3),  and  Equation  (4)  expresses 
conservation  of  air-mass  in  the  bag.  Collectively  they  are,  when  combined 
with  Equations  (1)  and  (5)  to  (11),  a  system  of  three  nonlinear,  ordinary 
differential  equations  for  the  three  functions  y,  V  and  p.  The  initial 
conditions  are 

at  t  =  0,  y=H,  V=  -Vo  and  p  ~  p  . 

a 

The  equations  can  be  put  in  dimensionless  form  by  defining 
xi  =  y/H,  x 2  =  V/Vo,  x 3  =  p/pa,  t  =  TH/Vo 
ai  =  paRb//W '  a2  =  gH/Vo2,  as  =  2ypa/[  (y-1  )pgVo 2  ] 

=  Pc/Pa 

a5  -  oaAcVc2CD/(2W) ,  ae  =  Jo 
n  =  p/pa  Q  =  q/Vo  0>  =  Ay/Ag, 


-3- 
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n 


which  leads  to  the  system  of  equations 


dx 

l/dT 

=  X  2 

(12) 

dx 

2/dT 

=  a2Fg 

(13) 

dx 

a/dT 

=  —  (  X  2  X  3 

+ 

QCv$)/xi 

(14) 

where 

n 

=  X  3 

v,  a*.  = 

[1 

4  ( y- 

1  )/2 ) 

1-(1/Y) 

(15), (16) 

Fg 

=  -1 

+  a5X22 

+ 

ai(n- 

1 

1) 

(17) 

Q 

=  [CX6X22  +  a 

30 

V 

(18) 

0 

=  (n/x3)-i 

n  i  a» 

(19) 

=  (Y 

-1 ) { n/ou 

1- 

(1/Y  ) 

n  i  cu 

(20) 

and  the 

initial 

conditions 

are 

XI 

=  1, 

X  2  = 

-1 

,  X  3 

=  1 

at  T  -  0 

(21) 

Typical  values  for  the  physical  constants  in  the  analysis  are  listed  in 
Table  1  and  the  corresponding  dimensionless  parameters  ctj  (J  =  are 
given  in  Table  2. 


TABLE  1 


Typical  Values  for  Physical  Constants 


H 

V 

A 

W 

P 

P 

g 

A 

A 


o 

B 

a 

a 

c 


v 


initial  height  of  platform 
initial  velocity  of  platform 
cross  sectional  area  of  bag 
weight  of  platform  and  load 
air  pressure 
mass  density  of  air 
acceleration  of  gravity 
drag  area  of  canopy 
cross-sectional  area  of  vent 


3  ft 
30  ft/s 
10  ft2 
1000  lbs 
2117  lb/ft2 
.002  lb  s2/ft“ 
32.2  ft/s 
1000  ft2 
1  ft2 


TABLE  2 


Typical  Dimensionless  Parameter  Values 


ai 

=  21.17 

a2 

=  .1073 

01  3 

=  8233 

a4 

=  1.893 

a5 

=  .  9C 

n 

a6 

-  90  S0 

If  the  airbag  has  any  effect,  we  may  assume  that 
j  x i |  •  1  and  | x* |  <  3 

tr.roughout  most  of  the  motion  because  <  1,  as  <<  ai,  and  it  is  plausible  to 
neglect  the  term  in  (17)  that  contains  as,  although  some  inaccuracy  can  result 
near  T  =  0.  A  similar  argument  causes  us  to  neglect  the  term  in  (18)  that 
involves  a& ,  although  the  accuracy  may  be  questionable  near  T  =  0,  especially 
since  the  value  of  Jo  is  not  well-specified  beyond  saying  that  Jo  =  0(1). 


The  remainder  of  this  paper  will  be  based  on  the  equations  (12)  to  (16), 
the  modified  equations  (17)  and  (18), 

Fa  =  -1  +  ai(n  -  1)  (22) 

1 

<2  =  (a 30) 5  (23) 

and  ( 19 )  to  ( 21 ) . 

With  these  approximations,  the  equations  involve  three  dimensionless 
parameters,  ai,  a:  and  03.  We  can  see  at  a  glance  that  model  tests  of  the 
platform-airbag  system  will  encounter  scaling  problems  unless  the  model  tests 
are  carried  out  in  a  suitable,  artificial  atmosphere.  For,  if  a2  is  to  be  the 
same  for  model  and  prototype,  we  must  have  VH~^  the  same,  but  as  implies  that 
V  must  be  the  same  for  model  and  prototype  if  both  are  tested  in  the  same 
medium. 

The  energy  possessed  by  the  platform  can  be  found  by  observing  that  the 
potential  energy  is  Wy  and  the  kinetic  energy  is  WV2/(2g).  Then 

C  =  (total  energy )/( initial  energy) 

=  ( 2a2X 1  +  X2  2 )/ ( 2a2  +  1)  (24) 

The  purpose  of  the  airbag  is  to  decelerate  the  platform  so  its  vertical 
velocity  is  zero  when  the  platform  strikes  the  ground,  i.e.,  xi  =  X2  =  0  at 
same  time,  T^.  The  airbag  has,  therefore,  to  completely  dissipate  the  initial 
energy . 


3.  SOLUTIONS  FOR  SPECIAL  CASES 

In  general  the  system  of  equations  is  nonlinear  and  cannot  be  solved 
exactly  in  closed  form.  However,  in  certain  special  cases  exact  or 
approximate  solutions  can  be  found,  and  it  is  convenient  to  describe  these 
here . 

It  is  instructive  first  to  think  qualitatively  about  the  behavior  of  the 
system.  For  this  purpose,  hodograph  plots  of  three  more  or  less  typical  cases 
are  shown  in  Figure  1.  All  start,  of  course,  from  the  initial  point,  xi  =  1, 
X2  =  -1.  The  path  Po  is  the  trajectory  followed  in  an  ideal  case  where  the 
system  reaches  the  point  xi  =  X2  =  0  at  time,  Tj,  possibly  under  the  action  of 
some  type  of  control  system.  A  case  where  the  bag  is  vented  too  freely,  and 
the  platform  crashes  into  the  ground,  i.e.,  X2  <  0  when  xi  =  0,  exhibits  a 
trajectory  like  P_.  Contrarily,  P+  shows  a  path  when  the  bag  is  not  vented 
enough,  and  the  platform  bounces  off  the  bag,  i.e.,  X2  =  0  when  xi  >  0. 


Figure  1:  Hodographs  for  platform-landings  in  three  typical  Cases: 

Optima]  Venting,  P0 ,  Overventing,  P  ,  Underventing,  P+. 
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The  first  special  case  is  that  in  which  <t>  =  0,  i.e.,  the  vent  is  closed. 

An  exact  solution  is  obtainable  in  implicit  form  in  this  case, 
f(xi)  =  {1  +  2a; ( 1+ai ) ( 1-xi )  +  2aia2(y-l)_1  (l-xi1_Y)} 

xi'1  h  =  xfY  (29) 


X  2  = 

-f’(Xl)  X  3 

1 

T  ( x  1 ! 

1  =  J  f”*(u)du 

This  solution  resembles  that  of  the  trajectory  P+  in  Figure  1  if  it  is 
followed  long  enough.  The  above  solution  loses  validity  when  the  function 
f(xi)  =  0,  i.e.,  when  X2  =  0.  The  values  of  xi  and  T  at  which  this  occurs 
depend  only  on  ai,  a2  and  y ,  not  on  as. 


An  important  special  case  is  obtained  if  we  demand  that  the  g-force,  F  , 
be  held  constant  from  some  time  (say  To  >  0)  until  the  final  time  T^.  For  the 
moment  we  ignore  what  happens  for  0  £  T  <  To. 


If 


xi(To)  =  bi  >  0,  X2(To)  =  b2  <  0,  (30) 

it  is  easily  verified  that  the  following  functions  not  only  satisfy  the  thre  ; 
equations  of  motion  but  also  the  end  conditions,  xi  =  X2  =  0  at  T  =  T^: 

Tf  =  To  -  2bi/b2,  Fg  =  b22/(2b1a2)  (31) 

(32) 

(33) 
(34), (35) 

<t>  =  -X2X3/Q  (36) 

and  Q  is  given  by  (18)  and  (19)  or  (20).  This  solution  has  constant  values 


xi  =  01  +  b 2 ( T-To )  +  ia2F  (T-To ) 2 

y 

X2  =  b2  +  a2Fg(T-To) 


X3  =  {1  +  ( 1+Fg )/ai } 


l/y 


n  =  X3 1 


for  F  ,  x;,  n  and  Q,  quadratic  time-dependence  of  xi,  and  linear  time 
y 

dependence  for  X2  and  $ .  Also  <j>  =  0  at  T  =  T^.  In  effect.  Equation  (36) 
describes  the  time-dependence  of  vent-opening  that  is  needed  in  order  to 
obtain  this  motion. 


It  is  clear  that,  if  this  motion  can  be  obtained  in  practice,  it  is  an 


optimum  solution  to  the  problem,  at  least  for  To  S  T  £  T^. 
section  we  shall  discuss  some  aspects  of  this  question. 


In  the  next 


4 .  CONTROL  SYSTEMS 


All  of  the  control  systems  studied  in  this  report  use  <J>  (the  vent-area 
ratio)  as  the  control  variable.  For  simplicity  in  studying  the  action  of  the 
control,  we  shall  also  assume  that  Cv  =  1  in  Equation  (14),  a  condition  that 
causes  some  loss  of  generality  if  Cv  depends  on  the  state  of  the  system  but 
not  otherwise. 

The  control  systems  that  have  been  studied  or  used  in  the  past  are  as 
follows : 

(i)  The  constant  vent  opening, 

0  = 

is  probably  the  simplest  system.  It  involves  only  a  single  parameter,  4>c. 

(ii)  The  blow-off  patch  has  the  vent  remaining  closed  until  the  pressure 
first  attains  a  certain  value,  Hq#  at  which  time  the  patch  blows  off 
(instantaneously)  and  the  vent  area  jumps  to  its  fully  open  value,  <(>c,  and 
remains  at  that  value  for  the  rest  of  the  landing,  i.e., 

-0  =  0  T  <  Tb 

<P  =  *c  T  >  Tb  (37) 

n(TB)  =  nB 

Thus  control  depends  on  two  parameters,  Hg  and  <j>c. 

Another  control  system,  not  previously  examined  nor  used,  but  apparently 
worth  investigation,  is  implicit  in  the  solution  of  equations  (30)  to  (36). 
This  solution  is  optimal  for  To  ST  S  Tj,  but  the  problem  is  in  determining 
when  to  start  using  the  control  law  (36)  which  causes  the  system  to  follow 
this  solution.  This  suggests  that  we  use  a  two-stage  system,  in  which  the 
first  stage  controls  <p  so  that  Fg  is  brought  quickly  from  its  initial  value 
(Fg  =  -1)  to  the  value  given  by  (31).  When  Fg  attains  that  value,  the  second 
stage  begins,  in  which  <f>  simply  follows  the  control  law  (36).  This  requires 
that  in  the  first  stage  the  system  must  sense  Fg,  xi  and  X2  and  switch  when 
(31)  is  satisfied,  i.e., 

Fg  =  X2  2/(2x ia2 )  (38) 
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xi  a.iu  x 2  can  be  calculated  by  sensing  Fg  and  performing  numerical  integration 
of  (12)  and  (13),  so  that  only  Fg  really  has  to  be  acquired. 

However,  another  question  is  whether  this  control  system  can  respond 
quickly  enough  so  that  4>  will  be  set  "instantaneously"  (i.e.  in  negligible 
time)  to  the  value  demanded  by  (36).  Alternatively,  can  $  be  controlled  in 
the  first  stage  so  that  it  will  both  (a)  pressurize  the  bag  enough  to  cause 
(31)  to  be  satisfied  at  some  To  and  (b)  have  the  value  demanded  by  (36)  at 
that  To? 

It  is  not  clear  whether  these  difficulties  can  be  surmounted  in  practice 

without  making  the  system  distastefully  complicated.  Also,  when  this  system 

is  in  the  second  stage,  the  automatic  control  is  of  open  loop  type,  i.e.,  it 

makes  no  use  of  information  about  F  or  any  of  the  stated  variables.  The  same 

y 

is  true  of  the  constant  vent  opening.  The  blow-off  patch  is  affected  by  the 
system  state  only  to  the  extent  of  being  actuated  when  n  is  large  enough.  To 
some  degree,  therefore,  these  systems  all  share  the  usual  shortcomings  of  open 
loop  control,  in  particular  they  may  not  function  well  if  subject  to  unknown 
or  random  fluctuations  of  input. 

Accordingly,  an  entirely  different  control  system  was  investigated  in  this 

study.  Since  measurements  of  F  are  the  easiest  ones  to  obtain,  the  system 

y 

was  defined  by 


D. 

=  P2Fg  -  Pi  =  P2 ( Fg-r ) ,  r  =  P1/P2 

(39) 

d<J)/dT 

=  if  >  0  or  d  >  0 

(39) 

=  C  if  D,  s  0  and  <p  -  0 

(39) 

where  the  parameters  Pi  and  P2  have  to  be  chosen  so  that  the  system  is  brought 
to  xi  =  X2  =0.  This  system  is  conceptually  very  simple;  the  vent  is  opened 
or  closed  at  a  rate  proportional  to  F  -r. 


5.  NUMERICAL  ANALYSIS  OF  LANDINGS  WITH  CONTROLS 


In.  order  to  investigate  the  behavior  of  the  control  systems  described  in 
Section  4,  a  computer  program  was  written  that  carries  out  the  numerical 
integration  of  the  differential  equation  systems  (12)  to  (16)  and  (19)  to 
(23).  This  program  is  listed  in  Appendix  A.  It  consists  of  a  main  program, 
MAIN,  which  reads  the  physical  parameters,  forms  the  dimensionless  quantities 
oi,  ct.2,  03,  and 

a“  =  nc, 

and  invokes  the  IMSL  routine  DVERK  to  do  the  numerical  integration.  After 
completing  the  solution,  MAIN  writes  the  results  for  xi,  x 2,  X3 , 

X4  = 

and  F  in  file  7  and  the  derivatives,  dx-/dT  (j  =  1,...,4),  in  file  8. 
y  j 

The  routine  DVERK  uses  a  Range-Kutta  integration  scheme  which  requires 
that  the  user  furnish  a  subroutine,  FCN,  for  evaluating  each  dx^/dT,  given  the 
values  of  the  time  and  all  the  Xj.  The  control  system  is  modelled  by  a  set  of 
instructions  in  this  subroutine  which  define  xi,  =  <J>  or  dxi,/dT  =  d<fi/dT. 

We  adopted  the  parameters  given  in  Table  2  as  a  standard  condition  and  for 
this  condition  explored  the  behavior  of  the  three  control  systems  described 
earlier,  recalling  that  the  objectives  are  to  bring  the  system  from  the 
initial  point 

Xl  =  1,  X  2  =  -1,  X3  =  1 

to  the  final  point 

XX  =  X2  =  0 

and  to  do  so  as  smoothly  as  possible  in  the  sense  that  the  maximum  of  Fg 
during  the  motion  should  be  as  small  as  possible. 


The  simplest  control  system  is  that  with  constant  vent  opening, 

x-  M  =  $c, 

or 

dxn/dT  =  0  and  xi.(0)  =  <t>c . 

With  this  single  parameter  it  was  impossible  to  attain  the  final  point 
Xi  =  X2  =  0.  The  "best"  results  (best  in  the  sense  described  below)  were 


found  from  a  series  of  trials  to  be  at  x<.  =  #  =  0.019.  For  this  case  the 

behavior  of  xi  and  X2  is  shown  in  Figure  2  and  that  of  F  in  Figure  3. 

y 


Figure  2:  Dimensionless  Height,  X:,  and  Velocity,  X2,  as  functions  of 

dimensionless  time,  T,  for  the  standard  condition  and  constant 
vent  opening,  0  =  0.019. 


Figure  3:  G-force  or  dimensionless  acceleration,  Fg,  as  a  function  of 

dimensionless  time,  T,  for  the  standard  condition  and  constant 
vent  opening,  0  =  0.019. 


The  system  has  X2  -  0  when  1.0  <  T  <  1.5,  but  xi  =  .36  at  that  time. 

After  almost  coming  to  rest  the  system  gradually  resumes  its  descent  until  it 
strikes  the  ground  (i.e.,  xi  =  0)  with  velocity  X2  =  -.18  at  time  T  -  4.1. 
During  this  second  phase  of  the  landing,  the  pressure  of  the  air  in  the  bag  is 
not  great  enough  to  equilibrate  the  load,  so  Fg  is  slightly  negative.  If 
x<4  =  $  >  0.019,  i.e.,  the  vent  is  opened  wider  than  the  optimal  value,  the 

load  strikes  the  ground  earlier  and  with  higher  velocity.  For  example,  if 
<}>c  =  0.025,  the  load  lands  when  T  =  2.17  with  X2  *  -.25.  On  the  other  hand  if 
<PC  <  0.019,  the  vent  is  narrower  than  the  optimal  value,  the  load  bounces  off 
the  bag  and  accuracy  of  Equations  (18)  to  (20)  is  thereafter  doubtful. 
Moreover,  the  maximum  Fg  values  are  higher  in  this  case  than  the  others. 

For  example 

<()  =  .015  causes  max  F  -  16 

<PC  =  .019  causes  max  Fg  =  14 

=  .025  causes  max  F_  =  9 

To  summarize,  this  simple  control  cannot  steer  this  system  to  the  origin. 
At  best  it  will  land  this  system  with  velocity  X2  =  -.18  and  max  Fg  =  14. 

The  blow-off  patch  control  system  was  examined  next.  Control  now  depends 
on  two  parameters 

nB  =  blow-off  pressure  (in  atmospheres) 

<(>c  =  dimensionless  vent  area. 

A  number  of  cases  were  run  for  various  values  of  these  two  parameters. 

The  results  were  qualitatively  like  those  for  the  previous  constant  vent 
control,  in  that  the  system  could  not  be  steered  to  the  point  xi  =  X2  =  0  by 
the  control.  Instead  the  system  attained  X2  =  0  at  height  xi ,  (i.e.  it  paused 
at  height  xi)  and  then  eventually  attained  xi  =  0  at  a  velocity  X2,  just  as 
for  the  constant  vent  case.  The  "best"  results  are  listed  in  Table  3. 


-12- 


TABLE  3 


nB 

Optimal  Results 

for  Blow- 

Off  Patch 

Tf 

*c 

*i 

><2 

F 

g 

1.3 

.019 

.40 

-.18 

14 

4.4 

1.5 

.44 

-.19 

14 

4.6 

1.8 

.023 

.47 

-.22 

16 

4.4 

Xj  is  the  value  of 

x1  at  which  x2  =  0  (height  of 

pausing) . 

x2  is  the  value  of 

x2  at  which  x1  =  0  (velocity  at  landing) 

F  is  the  maximum  value  of  F 

. 

g  g 


The  feedback  control  system  (39)  displayed  behavior  quite  different  from 
the  other  two  controls.  Many  different  pairs  of  values  for  Pi  and  P2  were 
found  that  would  steer  the  system  to  xi  =  X2  =  0.  For  example.  Table  4  shows 
three  such  sets  of  values  and  the  principal  properties  of  the  trajectories 
that  they  produced.  Figure  4  shows  graphs  of  xj  and  x:  as  functions  of  T  for 
the  case  Pi  =  1.02,  P2  =  .2,  and  the  functions  X3,  F  and  <p  are  depicted  in 
Figures  5,  6  and  7.  The  results  were  qualitatively  similar  for  the  other 
cases  listed  in  Table  4. 


TABLE  4 


Pi 

Computed  Results 

for  Feedback 

Control 

Law 

Tf 

?2 

F 

g 

M 

d$/dT 

Sf 

.50 

0.1 

9.3 

4 

-.027 

4.8 

1.95 

1.02 

0.2 

8.1 

6 

-.024 

4.9 

1.93 

2.08 

0.4 

7.4 

8 

-.026 

5.1 

1.90 

4.16 

0.8 

6.4 

9 

-.027 

5.2 

1.90 

F  is  the  maximum  value  of  F  . 

M*3  is  the  number  of  local  maxima  of  F  in  0  £  T  =  T 
d<J>/dT  is  the  approximate  slope  of  $  niar  T  =  T 

is  the  approximate  constant  value  of  G^.  near  T  =  T^ 
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The  principal  feature  of  these  results  is  the  oscillation  in  F  ,  X3  and  <f>, 

y 


an  oscillation  which  is  mildly  discernible  also  in  the  plot  of  X2  but  not  of 
xi.  A  perturbation  analysis  of  the  dif f erential-eguation  system  is  done  in 
Appendix  B  to  show  the  origins  of  this  behavior.  However,  it  is  clear  that 


the  F  values  obtained  with  this  control  are  much  lower  than  with  either  of 
9 


the  other  controls.  A  practical  question  is  whether  the  control  system  can 
respond  quickly  enough  to  enforce  the  control  law  during  an  oscillation  of 
this  type. 


TIKE,  T 


Fiaure  4: 


Dimensionless  Height,  Xlf  and  Velocity,  X2  as  Functions  of 
Dimensionless  Time,  T,  for  the  Standard  Condition, 

Feedback  Control  with  P2  =  1.02,  P2  =  0.20  and  T  2  1.0. 


TWE,  7 

Figure  5a:  Dimensionless  Air  Density  in  Bag,  X3,  as  a  Function  of 
Dimensionless  Time,  T,  for  the  Standard  Condition, 
Feedback  Control  with  Pt  =  1.02,  P2  =  0.20  and  T  £  1.0 
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Figure  5b:  Dimensionless  Air  Density  in  Bag,  X3,  as  a  Function  of 
Dimensionless  Time,  T,  for  the  Standard  Condition, 
Feedback  Control  with  P.  =  1.02,  P-,  =  0.20  and  T  £  1.0 


Figure  6a:  G-force  as  a  Function  of  Dimensionless  Time,  T,  for  the 
Standard  Condition,  Feedback  Control  with  P2  =  1.02, 

P2  =  0.20  and  T  S  1.0. 
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Figure  6b:  G-force  as  a  Function  of  Dimensionless  Time,  T,  for  the 
Standard  Condition,  Feedback  Control  with  P2  =  1.02, 

P,  =  0.20  and  T  2  1.0. 


6.  DISCUSSION 


The  solution  obtained  in  Appendix  B  agrees  quite  well  with  that  obtained 
by  the  computer  program  of  Appendix  A.  Superficially,  the  solutions  for  X3 
and  Fg  have  the  form  of  damped  oscillations  about  constant  values  with  periods 
that  decrease  as  T  -*•  T^,  and  this  is  exactly  what  is  seen  in  Figures  5  and  6. 
It  is  readily  verified  that  the  asymptotic  values 
x 3  ~  C3  =  1.195 
Fg  “  Pi/P2  ^  5 

agree  well  with  these  Figures.  Similarly,  the  asymptotic  result 
d$/dT  ~  -CL2Fg°C3/qo  =  -.0260 

conforms  closely  to  the  estimates  in  Table  4,  which  were  obtained  graphically 
from  Figure  7. 

To  assess  the  quantitative  agreement,  we  refer  to  Table  5,  which  shows  the 
local  maxima  of  X3  and  A3  and  their  times  of  occurrence,  .  Two  comparisons 
are  relevant.  First,  if  the  solution  (B.26)  is  correct, 

\  (In  t.j  -  In  ti)  -  2 it 

Ti-1/Ti  =  e2ll/X  =  e'256  =  L292- 

and  we  see  from  Table  5  that  successive  ratios  of  agree  very  well  with  this 
estimate.  Second,  the  maxima  should  satisfy 
In  A3i  =  A*  +  ( v-3/2 )  In  xi 

and  so 

d(lnA3i)/dlnTi  =  v  -  (3/2)  =  1.40. 

The  points  of  Table  5  give  a  value  approximately  =  1.50,  hence  there  is  a 
small  discrepancy  in  this  comparison. 
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TABLE  5 


T 

i 

OSC  1  I  lei*  I1 

on 

Extremes  for 

Feedback 

Control  Law 

luT  , 

1 

Ti  '  Tf  ‘ 

T. 

1 

T  .  /T  . 

1  1  +  1 

^3i 

InA , 

.  2b5 

1.665 

1.296 

.099 

-2.31 

.510 

.  64  6 

1.285 

1.298 

.065 

-2.73 

.251 

.940 

.990 

1 .303 

.044 

-3.12 

-.010 

1  .  170 

.760 

1 .288 

.030 

-3.51 

-.274 

1  .  340 

.590 

1.297 

.020 

-3.91 

-.528 

1.475 

.455 

1.319 

1  .  585 

.345 

1.302 

1.665 

.265 

Two  possible  sources 

of 

error  in  the 

estimate 

of  Appendix  B 

are  these 

(i)  The  series  expansions  involved  in  obtaining  Equations  (B.10),  (B.ll) 
and  (B.12)  are  not  extremely  accurate. 

<ii)  The  estimate  (B.23)  is  only  a  moderately  accurate  approximation  to 
the  solution  of  (B.22). 

It  is  questionable  whether  the  effort  involved  in  improving  these 
approximations  is  worth  the  trouble. 

Several  other  comments  can  be  made  about  this  solution. 

(a)  The  differential  equation  (B.25)  has  a  singularity  at  T  =  0,  T  = 
that  casts  doubt  upon  the  correctness  of  the  limiting  behavior  of  the 
solution,  (B.26),  as  -*■  T^. 

(b)  The  fact  that  F  0  =  Pi/P?  furnishes  a  fairly  accurate  solution  of 

y 

(B.22)  suggests  that  the  aero-th  order  solutions  are  not  very  sensitive  to  the 
values  of  Pi  and  P2  separately,  but  only  to  the  ratio  of  Pi/Pj,  a  conclusion 
supported  by  the  numerical  results  in  Table  4.  However,  the  perturbation 
solution  (B.27)  depends  strongly  on  A,  which  depends  on  a  and  so  is  influenced 
by  P2  alone  (see  (B.20)).  As  P2  increases,  A  increases  and  the  solutions 


oscillate  more  rapidly,  a  result  that  is  confirmed  by  the  values  of  M  in 
Table  4. 


(c)  A  valid  question  about  this  control  system  is  whether  the  control 
mechanism  can  change  the  vent  opening  fast  enough  to  keep  up  with  the  changes 

in  F  .  We  have 

y 

j  d<f/dT  |  £  |  d<$  o  /dT  |  +  |  dAo  / dT  j 

and 


difo/dT  =  -.026.  Also,  from  (B.20) 

| dAo/dT |  =  |aAa |  £  a) A3 | 

and  from  the  computations.  Figure  5,  we  see  that  [ A 3  j  £  0.1. 
Hence 


|  d<J)/dT  |  £  .026  +  6.37  x  .01  =  .090. 

Thus 

dAv/dT  =  VdAqH'1  d<J>/dT  £  9  ft2/s 

This  rather  crude  estimate  implies  that  the  system  must  be  able  to  change  the 
vent  area  at  a  rate  of  .009  ft2/ms  in  order  to  control  the  air  bag  in  the 
manner  assumed  by  the  analysis.  It  is  not  known  whether  this  is  an  attainable 
rate  because  much  depends  on  the  shape  of  the  vent,  but  for  any  specified  vent 
geometry  (e.g.,  rectangular),  it  should  be  possible  to  decide  the  question. 


(d)  We  see  from  Table  4  that,  as  Pi  and  Pz  increase  but  remain  in  almost 
the  same  proportion,  the  maximum  F^  decreases  noticeably.  This  suggests  that 
the  most  effective  control  is  obtained  with  P1/P2  -  5  and  P2  as  large  as 
possible.  However,  as  mentioned  under  (b),  the  oscillation  becomes  more  rapid 
as  P2  increases.  When  P2  is  increased,  the  oscillation  eventually  becomes  so 
quick  that  the  control  system  cannot  keep  up  with  it,  and  the  present  analysis 
becomes  inaccurate.  A  more  perceptive  analysis,  in  which  the  effect  of 
control  system  response  is  modelled  both  theoretically  and  numerically,  would 
shed  valuable  light  on  the  practical  improvement  that  might  be  attained  with  a 
control  system  of  this  general  type. 


(e)  The  parameter,  X,  that  determines  the  frequency  of  the  perturbation 
oscillation,  is  rather  large,  eg.,  X  ~  24.9  >>  1  in  the  example  of  Appendix  B. 
We  see  from  (B.26)  that  the  size  of  A  depends  on  qo  (i.e.,  ultimately  03),  a, 
and  r,.  While  qo  -  24.7  is  fairly  large,  a/b  ~  23.8  is  almost  as  large,  and 
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the  largeness  of  \  results  about  equally  from  both.  In  fact  a  is  fairly  large 

because  ai  is  so,  and  b  is  rather  small  because  ct?  is.  Thus  all  three 

”1  - 

parameters  cu,  aj  and  a3  in  the  combination  cua2  aa2  have  substantial 
influence  on  A. 

Concerning  the  two  open-loop  control  systems,  little  further  comment  is 
needed,  except  to  remark  that  the  numerical  integration  subroutine,  DVERK, 
experienced  some  convergence  difficulty  with  the  discontinuity  that  occurs  for 
the  blow-off  patch.  The  results  for  these  computations  are  less  accurate  than 
for  the  other  control  systems  although  not  sufficiently  so  to  alter  the  main 


7.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  of  this  preliminary  study  of  airbag  control  systems  are  as 
follows : 

(i)  Neither  of  the  open  loop  systems  performed  in  a  wholly  satisfactory 
way  for  two  reasons.  They  did  not  bring  the  system  to  the  desired  end  point, 
and  the  maximum  G- force  was  very  high. 

(ii)  The  feedback  system  with  vent-rate  proportional  to  the  G-force  was 
better  on  both  counts  than  either  open-loop  system.  For  many  pairs  of  values 
of  Pi  and  P2  it  brought  the  system  to  the  end  point  and  did  so  with  a  maximum 
G-force  much  smaller  than  either  open  loop  control. 

(iii)  A  number  of  questions  remain  about  this  particular  control  system. 

It  would  be  desirable  to  eliminate  the  oscillation  or  at  least  reduce  its 

frequency  if  that  can  be  done  without  seriously  increasing  F  or  degrading  its 

y 

ability  to  attain  the  desired  end  point.  Also,  we  need  to  clarify  the 
behavior  near  T  =  Tf,  and  it  would  be  desirable  to  conduct  stability  studies 
of  this  control  system,  i.e.,  how  it  responds  to  either  deterministic  or 
random  errors  in  the  inputs  or  environment. 

(iv)  The  results  encourage  us  to  think  that  closed-loop  control  systems 
in  general  have  much  to  offer  in  improving  air-bag  performance.  Although  the 
closed-loop  system  studied  in  this  report  is  a  plausible  one  that  may  be 
realizable  in  practice,  there  are  many  other  possibilities.  For  example  the 
control  law 

d4>/dT  =  -Pt  +  P2F  -  P3<tPl* 
y 

may  also  be  realizable  and,  if  the  parameters  are  chosen  well,  have  fewer 
undesirable  side  effects  than  the  law  studied  here. 

(v)  Concerning  the  computer  program,  the  values  of  the  control 
parameters,  Pi  and  P2,  were  found  by  manual  trial-and-error  in  the  present 
study.  It  is  possible  to  include  in  the  program  a  subroutine  for  nonlinear 
optimization,  which  will  carry  out  this  process  "automatically".  For  example, 
the  nonlinear  least-squares  solver  NL2S0L  has  a  number  of  attractive  features 


for  a  study  of  this  kind.  However,  even  the  best  of  these  subroutines  can 
fail  if  the  optimal  solution  is  not  unique,  so  much  caution  is  required  in 
their  use. 


We  recommend,  therefore,  that  further  study  of  both  open  and  closed-loop 
control  systems  be  undertaken,  with  emphasis  on  the  latter.  In  particular, 
further  study  of  the  present  control  algorithm  is  justified,  as  sketched  under 
(ill)  above.  Otner  control  laws  ought  also  to  be  examined.  If  a  substantial 
investigation  is  undertaken,  the  computer  program  of  Appendix  A  should  first 
be  enhanced  by  inclusion  of  a  carefully  chosen  nonlinear  optimization 
subroutine . 
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APPENDIX  A:  Computer  Programs 


The  following  pages  contain  listings  with  comments  of  the  three  elements 
needed  to  carry  out  the  computations  described  in  Section  5.  These  elements 
are  presently  stored  in  the  file 
BAG  *  FF. 

on  the  UNIVAC  1106  at  the  U.S.  Army  Natick  Research,  Development  and 
Engineering  Center.  The  elements  are 
M,  the  main  program 

FNI,  the  subroutine,  FCN,  called  by  DVERK 
D,  a  typical  data  set. 

M  is  listed  below  followed  by  FNI  and  D.  The  data  in  D  is  that  which  produced 
the  output  for  the  constant  vent-opening  control  system  in  Figures  2  and  3. 


PROGRAM  FOR  AIRBAG  CONTROL 


This  program  calculates  the  motion  of  a  platform  during  a  landing 
cushioned  by  an  airbag  with  an  automatic  control  system.  It  is  based  on  the 
system  of  dimensionless  equations  in  the  report  "A  Preliminary  Study  of 
Control  Systems  for  Platform  Landings  Cushioned  by  Airbags"  by  E.W.  Ross. 

The  main  program,  given  below,  reads  in  the  physical  quantities,  converts 
them  to  dimensionless  parameters  and  then  calls  the  IMSL  subroutine  DVERK, 
which  does  numerical  integration  of  the  differential  equation  system,  using  a 
Runge-Kutta  procedure.  The  subroutine  DVERK  call  the  user-supplied  subroutine 
FCN,  which  calculates  the  derivatives,  given  the  function  values.  The 
instructions  which  define  the  control  system  are  embedded  in  this  subroutine. 


DEFINITIONS  OF  PRINCIPAL  QUANTITIES: 


PA, RHOA,GM 
AB 
G 
W 


VO 

H 

N 

TOL 

NTS 

TI,  TF 

YI 

LAM 


CP 

Y 


DRV 

GF 


pressure,  mass-density  and  gamma  for  air 

cross-sectional  area  of  airbag 

acceleration  of  gravity 

weight  of  platform,  load  and  bag 

velocity  of  platform  at  first  contact 

height  of  platform  at  first  contact 

number  of  variables  in  vector  X  (usually  4) 

accuracy  threshold  for  DVERK 

number  of  time  steps  in  the  integration 

initial  and  final  time 

array  of  initial  values  of  the  variables  Y 

array  of  integer  variables  for  information  recording 
and  imposing  certain  conditions. 

array  of  control  and  information  paramaters 
array  of  main  variables,  as  follows: 

Y ( I )  =  dimensionless  height 
Y(2)  =  dimensionless  velocity 
Y(3)  =  dimensionless  density 

Y(4)  =  dimensionless  vent  area  (control  variable) 

array  of  derivatives  of  Y 

G-force  (i.e.,  acceleration  in  G-units) 


The  calculated  values  of  the  Y's  and  GF  are  written  to  file  no.  7  and 
DRV's  are  written  to  file  no.  8.  A  typical  set  of  input  data  for  use  with 
this  program  is  in  the  element  D  in  this  file.  The  program  graph  in  this  file 
can  be  used  to  cause  Tektronix  plotting  of  the  results  stored  in  file  7. 


-28- 


COMPILER  ( DIAG=3 ) 

EXTERNAL  FCN 

COMMON  /A/  GM , AL ( 4 ) , LAM ( 5 ) , GF , CP ( 5 ) DRV ( 4 ) 

REAL  YI(4) ,  Y(  4  ) , C ( 24 ) , WKV (10,10) 

****  READ  IN  PHYSICAL  AND  OTHER  DATA  ***** 

READ  (5,30)  PA , RHOA , GM 

READ  (5,30)  AB,G,W,UO,H 

READ  (5,30)  N,TOL,NTS 

READ  (5,30)  TI , TF, ( YI ( J) , J=1 ,4 ) 

READ  (5,30)  (LAM(K) ,K-1,5) 

READ  (5,30)  ( CP( J) , J=1 , 5 ) 

****  CALCULATE  THE  ARRAY  AL  (ALPHA)  OF  DIMENSIONLESS  PARAMETERS 
AL( 1 )=PA*AB/W 
AL(2)=G*H/UO/UO 

AL ( 3 )  =  2  *GM*  PA/ ( GM- 1 ) /RHOA/UO/UO 
AL( 4 ) = ( 1+GM- 1 )/2 ) ** ( GM/ ( GM-1 ) ) 

WRITE  (6,40)  ( AL( K) , K=1 ,4 ) 

IND=  1 
T=TI 

DO  10  J=1,N 

10  Y( J)=YI(J) 

DT=( TF-T) /NTS 

****  MAIN  LOOP  FOR  NUMERICAL  INTERATION 
DO  20  1=1, NTS 

TS=TI=I*DT 

CALL  DVERK  ( N ,FCN ,T , Y,TS,TOL, IND , C, 10 , WKV , IER) 

WRITE  (7,40)  T,  ( Y ( J ) , J= 1 , 4 ) , GF , CP ( 4 ) 

WRITE  (8,50)  T,  (DRV( J) ,J=1,4) 

20  CONTINUE 

****  END  OF  MAIN  LOOP 

30  FORMAT  ( ) 

40  FORMAT  (7E9.4) 

50  FORMAT  (F6.4,7E9.4) 

END 


SUBROUTINE  FCN 


This  subroutine  is  called  by  DVERK  to  evaluate  the  derivatives  (DY)  of  the 
variables,  given  their  values  and  those  of  the  paramameters .  The  control 
instructions,  starting  at  line  23,  are  for  the  particular  feedback  control 
system  studied  in  the  report,  which  involves  the  derivative  of  the  vent  area, 
Y(4).  For  a  different  control  system,  a  different  set  of  instructions  may 
have  to  be  written  and  perhaps  even  inserted  at  a  different  point  in  the 
program. 

COMPILER  ( DIAG=3 ) 

SUBROUTINE  FCN  ( N , T , Y , DY ) 

REAL  Y( N ) , DY ( N ) 

COMMON  /A/  GM , AL ( 4 ) , LAM ( 5 ) , GF , CP ( 5 ) , DRV ( 4 ) 

DY  fl)=Y(2) 

IF  ( Y ( 3 ) . LT .  0.0)  WRITE  (6,20)  T, (Y( J) ,J=1,4) 

PR=Y ( 3 ) **GM 
FD=LAM( 1 ) *Y( 2 ) **2 
BR=AL( 1 ) * ( PR- 1 . ) 

GF=- I . +FD+BR 
DY( 2 ) =AL( 2 ) *GF 
LAM( 5 ) =0 

AA=AL( 3 ) * ( PR/Y( 3 ) - 1 .  ) 

AP=AL ( 3 ) * ( GM- 1 ) /2 

IF  ( PR. GTAL( 4 ) )  AA=AP* ( PR/AL( 4 ) )**( 1- l/GM ) 

IF  (PR.GTAL(4) )  LAM(5)=1 

IF  (AA.LT.  5.E-8)  AA=0 . 

Q=S0RT( LAM( 2 ) *Y( 2 ) **2+AA) 

DY(3)=-(Y(2)*Y(3)+Q*Y(4))Y(1) 

*****  THE  FOLLOWING  INSTRUCTIONS  EXTERT  CONTROL 
D4=CP( 2 ) *GF-CP( 1 ) 

DY( 4 ) =D4 

IF  (  (D4.LT.0.0) .AND. (Y(4) .LE.0.0)  )  DY(4)=0.0 

DO  10  J=l,4 
10  DRV( J)=DY( J) 

RETURN 

20  FORMAT  (5E9.4) 

END 

*****  DATA  SET  D  ***** 

2117. .  .002.1.4 

10. . 32.2.1000. .30. .3. 

4,1. E-4 , 40 

0.0, 4.0,  1.0,-1.0,1.0,0.019 
0.0. 3,0, 1 

0.00,0.00,0. ,0. ,0. 
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APPENDIX  B:  Perturbation  Analysis 


This  appendix  presents  a  perturbation  solution  of  the  motion  equations 
with  feedback  control.  This  analysis  is  motivated  primarily  by  the  observed 
results  of  the  computations  and  secondarily  by  the  exact,  optimal  solution  in 
Section  3  for  Fg  =  constant. 


The  basic  set  of  equations  is  (12)  to  (  15 ) ,  ( 19 ) , ( 22 ) ,  ( 23 )  with  Cv  =  1  and 


the  control  law  (39).  The  equations  are  then 

dx1/dT  =  x  2  ;  dx2/dT  =  ct2Fg 

d(xixj)  =  -Qd  ;  d(p/ dT  =  P2F  -  Pi 

y 


Q  =  aa2(X3'!r”'1'-l)  5 
F  =  -1  +  ai(X3>-1-l). 

y 


(B.1,2) 

(B.3,4) 

(B.5) 

(B.6) 


mi 


nil 

o 


We  are  assuming  that  the  vent  is  open  and  the  flow  through  it  is  subsonic, 
hence  this  approximate  solution  is  not  valid  initially,  when  the  vent  is 
closed.  Figures  5  and  6  show  that  for  T  2  0.3  the  variables  xa  and  Fg 
oscillate  with  decreasing  amplitude  about  constant  values,  and  this  is  the 
beha"ior  that  we  seek  to  explain  from  the  above  equations. 


We  assume  that 

xj  =  xj°  + 
$1  =  <fo  + 


j  =  1,2,3 
F„  =  F  °  + 


(B .  7 ) 
(B.8,9) 


The  quantities  A-,  A,,  Ap  are  assumed  to  be  perturbations  of  the  zero-th  order 
0  •*  v 

quantities  x •  ,  <po,  F  °,  respectively,  with  A-  <<  xH 0 .  Series  expansions  of 

j  y  j  j 

(B.5)  and  (B.6)  lead  to 

Q  -  qo(l  +  qiA 3 )  (B.10) 

qo  =  a3M(x3°)T_1-l]^,  q  1  =  i(y  - 1 ) (x3 0 )y~2/ [ (x3 0 )y~1-l ]  (B.ll) 


0 \Y~2, 


1  0  1 


Fg°  =  ai[(x3V-1-l]-l  , 


0  y-1 

Av  =  a. iy  ( x 3  )  ’  A3  . 


(B. 10) 
(B.ll) 
( B . 12 ) 


m 


Equation  (B.12)  shows  that  F  0  is  constant  if  x3°  is  so.  Since  we 

y 

expect  to  obtain  a  solution  such  that  x3  and  Fg  oscillate  about  constant 

values,  we  assume 
„  0 


and  so 


=  C3  =  constant 


Fg°  =  ai(C3Y_1-l )-l. 


( B . 13  ) 


( B . 14  ) 


Then  (B.l)  and  (B.2)  imply 

xi°  =  b( T-Tf ) 2  ,  x 2 0  =  2b( T-Tf ) 
b  =  *a2Fg° 

Equations  (B.3)  and  (B.4)  become,  using  (B.12) 

d[(xi°  +  A i ) ( C3  +  A  3  )  ]  /  dT  =  -q°  ( 1  +  qiAsH^o+A^) 
d(4>°  +  A^/dT  =  P2[Fg°  *an  C3y"1A3]-Pi. 

The  zero-order  terms  lead  to 
C3X20  =  -qo<J>o 
d<t>o/dT  =  P2Fg°  -  Pi 
and  the  first  order  terms  to 

d(xi °A3 )/dT  =  -C3  dAi/dT  -  qofA^  +  qi<fioA3) 
dA^/dT  =  aA3  ,  a  =  P2011Y  C3Y_1. 

Equations  (B.17)  and  (B.15)  imply 

<t>°  =  -2bC3(T-Tf  )/qo  (B.  21 ) 

and  ( B . 18 )  becomes 

P2Fg°  -  Pi  +  2bC3/qo  =  0.  (B.22) 

Since  Fg°  and  qo  all  depend  on  C3 ,  this  is  a  transcendental  equation  for 
C3  and  has  to  be  solved  by  trial  and  error. 

However,  we  see  from  (B.ll)  that  qo  involves  the  parameter  a3  which  is 
very  large  for  the  present  set  of  parameters,  see  Table  2.  This  implies  that 
qo  >>  1  and  suggests  that  we  attempt  to  solve  (B.22)  by  neglecting  the  last 
term,  which  leads  to 

Fg°  =  PVP2.  (B.23) 

For  all  the  values  of  Pi  and  P2  in  Table  4,  i.e.  those  values  for  which 


(B.17) 
(B. 18) 

( B . 19 ) 
( B . 20 ) 


(B.15) 
( B . 16 ) 


Xl=  x  2 


0  at  T  =  Tf,  we  have 


Wv  w- J 


This  implies  via  (B.14),  (B.ll)  and  (B.16) 

X3J  =  C3  =  li  +  (  1  +  Fg°  )/oi  ]1/A(  =  1.195 
q b  -  24.7,  ui  ^  2.43  and  b  -  .268, 

and  we  car.  verify  that  the  last  term  in  (B.22)  does  not  greatly  affect  the 
estimate  (B.23). 


These  results  are  now  used  to  solve  (B.19)  and  (B.20).  We  can  eliminate 
A^  by  differentiating  (B.19)  and  substituting  (B.20),  obtaining 

d2  (Xi°A3  )/dT2  +■  qoq i(J> odA 3 / dT  +  C3  c^/dT  +  a  qoA3  +  qoqiAsdcJio/dT  =  0. 


With  the  aid  of  (B.2),  (B.7)  and  (B.12)  we  find 
C3dAj/dT  =  PA3,  p  =  aia2yC3^. 

From  (B.21) 

qaqid<f>o/dT  =  -2bv,  v  =  C3qi. 

Finally,  if  we  define 
u  =  x  1 0  A  3 

and  use  (B.15)  we  get 

d2u  -  2  vx_1du  +  (A2+2v)t~2u  =  0 


( B . 24  ) 


(B.25) 


where 


t  =  Tf  -  T  i  0 
A2  =  (ago  +  p)/b. 


( B . 26 ) 


The  numerical  values  of  these  quantities  are  in  this  case 
v  =  2.90,  p  =  4.08 

and,  for  the  case  where  P2  =  .2,  line  2  of  Table  2, 
a  =  6.37,  A  =  24.5. 


A  general  solution  of  this  equation  is 
u  =  At^v+^  cos(3  In  t  -  S) 

P  =  {A2  +  2v  -  (l+2v)2/4}2  =  A  +  0( A-1 ) 
where  A  and  S  are  the  arbitrary  amplitude  and  phase. 


The  solution  for  A3  is  \ound  from  (B.24)  and  (B.15) 
A3  =  AW*"*3/2)!  cos(  P  In  t  -  S). 

Also  Ap  and  dA^./dT  can  be  found  from  (B.12)  and  (B.20). 


(B.27) 


Kb  <1 


APPENDIX  C:  List  of  Symbols 


i 


A, A1 
AB,Av 


a 


b 


C 


D 


v 

C3 


D 

F, 

g 

H 


Jo 

M 


P 

P  1  ,  P  2 

Pa 

Pc 

Q 

q 

qo  ,qi 
rb 

r 

S 

t 

T 

Tf 

To 

u 

V 

Vo 

W 


Xl  ,X2  ,X3 


X  I  0  ,  X  2  °  ,  X  3  0 


Arbitrary  constant  in  perturbation  solution 
Cross-sect  ion  areas  of  airbag  and  vent 
Drag  area  of  canopy 

Parameter  in  perturbation  analysis,  see  (B.20) 

Parameter  in  perturbation  analysis,  see  (B.16) 

Parachute  drag  coefficient 
Vent  flow  coefficient 

Constant  density  value  in  perturbation  analysis,  see  (B.13) 
Canopy  drag,  see  Equation  (3) 

G-force  or  dimensionless  acceleration,  see  Equation  (17) 
Acceleration  of  gravity 
Height  of  airbag 

Constant  in  air-flow  definition,  see  Equation  (8) 

Number  of  oscillations  (local  maxima)  of  Fg  in  Table  4 
Pressure  of  air  in  airbag 

Constants  in  feedback  control  system,  see  Equation  (39) 

Standard  atmospheric  pressure 

Critical  (sonic)  pressure  in  vent 

Dimensionless  air-speed  in  vent 

Air-speed  in  vent 

Constants  in  perturbation  of  Q,  see  (B.ll) 

Reaction  (lift)  of  the  airbag  on  platform,  Equation  (3) 
P1/P2 

Variable  in  the  vent  air-speed,  Equation  (8) 

Time 

Dimensionless  time 

Time  at  which  platform  strikes  ground 
Initial  time  for  optimal  solution 
Variable  in  perturbation,  see  (B.24) 

Platform  velocity  of  descent 

Initial  platform  velocity 

Weight  of  platform  and  load 

Dimensionless  height,  velocity  and  density 

Zero-order  perturbations  in  xi,X2,X3,  see  (B.7) 
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( 


A j , A  3 , Ap, 


Height  of  platform  above  ground  during  landing 
Dimensionless  parameters 

Constant  in  perturbation  solution,  see  (B.27) 

Ratio  of  specific  heats  of  air 
Perturbations  of  solution,  see  (B.8),  (B.9) 

Dimensionless  pressure  of  air  in  bag 
Pressure  at  which  blow-off  patch  is  activated 
Critical,  sonic  pressure  in  vent 
Energy  of  platform 

Large  parameter  in  perturbation  oscillation,  see  (B.26) 
Constant  in  perturbation  analysis 
Constant  in  perturbation  analysis 
Mass  density  of  air  in  bag 

Mass  density  of  air  at  standard  atmospheric  conditions 
Dimensionless  variable  in  vent  flow,  see  Equations  (18),  (19) 
Tf-T,  see  (B.26) 

Dimensionless  vent  area,  control  variable 
Constant  vent  opening 

Zero-th  order  perturbation  in  vent  opening,  see  (B.8) 


fyyt 


